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Abstract. Internal rotation is considered to play a major role in determining the structure 
and dynamics of some globular clusters. We present a dynamical analysis of the photometry 
and three-dimensional kinematics of 47 Tuc and a> Cen, by means of a new family of self- 
consistent axisymmetric rotating models. The combined use of line-of-sight velocities and 
proper motions allows us to obtain a global description of the internal dynamical structure of 
the objects together with an estimate of their dynamical distances. The well-relaxed cluster 
47 Tuc is very well interpreted by our dynamical models; in particular, internal rotation 
is found to explain the observed morphology. For the partially relaxed cluster ai Cen, the 
selected model provides a good representation of its complex three-dimensional kinematics, 
in general qualitative agreement with the observed anisotropy profile, which is characterized 
by tangential anisotropy in the outer parts; discrepancies are found between the observed 
and the expected ellipticity profile and are ascribed to the presence of a high degree of radial 
anisotropy in the intermediate region and to its interplay with rotation. 
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1. Introduction 

Globular clusters (GCs) have long been con- 
sidered simple quasi-relaxed nonrotating stel- 
lar systems, characterized by spherical symme- 
try and isotropy in velocity space. Spherical 
nonrotating isotropic models (e.g., the King 
1966 and the spherical Wilson 1975 models) 
have indeed provided an adequate zero-order 
description of their internal dynamics. 

However, deviations from sphericity are 
observed and can now be measured in quan- 



titative detail (see Geyer et al. 1983, White & 
Shawl 1987, and Chen & Chen 2010). In addi- 
tion, significant internal rotation has been de- 
tected in a growing number of GCs from line- 
of-sight velocity measurements (Bellazzini et 
al. 2012) and, in a few cases, from proper mo- 
tion measurements (e.g., van Leeuwen et al. 
2000; Anderson & King 2003). Detailed three- 
dimensional kinematics is therefore available 
for selected GCs, in particular thanks to proper 
motion data sets from the Hubble Space 
Telescope (e.g., Anderson & van der Marel 
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2010). This progress calls for the development 
of a more realistic dynamical modeling frame- 
work, in which rotation and deviations from 
sphericity are taken into consideration. 

Interest in the study of internal rotation de- 
rives from a number of open problems: (1) in- 
ternal rotation, together with external tides, and 
pressure anisotropy, is considered to be one 
of the main physical factors responsible for 
the observed flattening of GCs (van den Bergh 
2008); (2) the presence of global angular mo- 
mentum is expected to change the long-term 
dynamical evolution of stellar systems (Fiestas 
et al. 2006); (3) difi'erential rotation can co- 
operate with pressure anisotropy to produce 
nontrivial gradients in the kinematic profiles 
(Varri & Bertin 2012). Unfortunately, in only 
few cases has internal rotation been studied by 
a quantitative application of nonspherical ro- 
tating dynamical models. 

Here we present, as a short preliminary re- 
port on an extensive study presented in a sepa- 
rate article (Bianchini et al. 2013, paper sub- 
mitted), the application of a family of self- 
consistent rotating models recently constructed 
(Varri & Bertin 2012) to two Galactic GCs, 
47 Tuc and cj Cen. The dynamical models 
are compared to the photometric and three- 
dimensional kinematic data. Furthermore, by 
taking into consideration the inclination angle 
of the rotation axis of the stellar systems, we 
perform a detailed analysis of the morphol- 
ogy of the clusters, thus testing its connection 
with rotation. Finally, since the two clusters 
are in difi'erent relaxation states (47 Tuc with 
logTc < 8 falls in the class of well relaxed 
clusters, whereas lo Cen with log Tc > 9 should 
be considered as only partially relaxed; here Tc 
indicates the core relaxation time expressed in 
years, Zocchi et al. 2012) we can test whether 
intemal rotation plays dilferent roles in sys- 
tems under difl'erent relaxation conditions. 



defined by a distribution function dependent 
on the integral of the motion I(E, J^) = E - 
{LoJy)l(\ + bJ^), such that the rotation is differ- 
ential, that is, approximately rigid in the cen- 
ter and declining and eventually vanishing in 
the outer parts. The velocity dispersion ten- 
sor is characterized by isotropy in the cen- 
tral region, weak radial anisotropy in the in- 
termediate regions, and tangential anisotropy 
in the outer parts (this behavior results from 
the requirement of self-consistency and from 
the adopted truncation prescription in phase 
space). Three dimensionless parameters define 
the models: the dimensionless depth of the po- 
tential well *P, the rotation- strength parame- 
ter, and the parameter h defining the shape of 
the rotation profile. Furthermore, 5 additional 
quantities must be specified: 2 physical scales 
(e.g., the radial scale and the velocity scale), 
the mass-to-light ratio (to convert mass density 
profiles into surface brightness profiles), the in- 
chnation angle of the rotation axis with respect 
to the Une of sight, and the distance to the clus- 
ter (to convert proper motions into km s"'). 

Models and observations are matched 
through a kinematic priority approach in the 
following steps. First, we determine the dimen- 
sionless parameters that characterize the inter- 
nal structure of the models by following natural 
criteria suggested by the observed kinematics, 
in particular by the characteristics of difl^eren- 
tial rotation. Then, we set the physical scales 
by means of standard fits to the surface bright- 
ness profile and line-of-sight dispersion and ro- 
tation profiles. An additional fit to the proper 
motion dispersion profiles gives us an estimate 
for the distance of the cluster Finally, once all 
the scales have been fixed, other observable 
quantities follow as predictions; in particular, 
the selected model makes a well defined pre- 
diction on the ellipticity profile. 



2. Dynamical models and fitting 
procedure 

The family of self-consistent axisymmetric 
models used in our analysis was specifically 
designed to describe quasi-relaxed stellar sys- 
tems (Varri & Bertin 2012). These models are 



3. 47 Tuc: a well-relaxed globular 
cluster 

The kinematic data sample for 47 Tuc con- 
sists of 2476 line-of-sight-velocities covering 
the entire extent of the cluster from Gebhardt 
et al. (1995) and Lane et al. (201 1) and 12 974 
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HST proper motions from McLaughlin et al. 
(2006) limited to the central 4 core radii. 

The photometric and kinematic profiles are 
well reproduced by the model, both along the 
line of sight and on the plane of the sky, con- 
firming isotropy in velocity space in the cen- 
tral region. In particular, the rotation profile is 
well matched throughout the radial extent of 
the cluster and the observed ellipticity profile 
follows the prediction of the model with sur- 
prising accuracy (Fig. 1). Given the fact that 
the selected model is associated with an ellip- 
ticity profile that is the morphological counter- 
part to the presence of rotation, we conclude 
that the observed deviations from sphericity 
originate from the presence of intemal rotation. 

Finally, we estimate a dynamical distance 
to 47 Tuc of d=4.15 ±0.07kpc, consistent with 
previous dynamical estimates, and slightly 
smaller than the distance derived from photo- 
metric methods (McLaughlin et al. 2006). 



4. oj Cen: a partially relaxed globular 
cluster 

The kinematic data sample for co Cen consists 
of 1868 line-of-sight- velocities from Reijns et 
al. (2006) and Pancino et al. (2007), reach- 
ing a radial extent of approximately half trun- 
cation radius, 2740 ground-based proper mo- 
tions from van Leeuwen et al. (2000), and 
72 970 HST proper motion measurements from 
Anderson & van der Marel (2010), with radial 
extent of half truncation radius. 

The surface brightness and the line-of-sight 
kinematic profiles are in satisfactory agree- 
ment with the model and the predicted average 
ellipticity is consistent with the one observed. 
However, discrepancies are noted for the ellip- 
ticity profile (Fig. 2). We argue that the sig- 
nificant ofl'set between observed and predicted 
ellipticity profile confirms the complex nature 
of OJ Cen, resulting from its condition of par- 
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Fig. 1. Line-of-sight rotation profile measured 
along the major axis (top panel) and ellipticity pro- 
file (bottom panel) for 47 Tuc. Solid lines represent 

the model profiles, open circles the observational 
kinematic data points, and black dots the observed 
ellipticities from White & Shawl (1987). 
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<o Cen - Ellipticity profile 




Fig. 2. Line-of-sight rotation profile measured 
along the major axis (top panel) and ellipticity pro- 
file (bottom panel) for u) Cen. Solid lines represent 
the model profiles, the open circles the observed el- 
lipticities from Anderson & van der Marel (2010), 
the black dots those from Geyer et al. (1983). 
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CO Cen - Anisotropy profile 
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Fig. 3. Anisotropy profile of co Cen, defined as the ratio of the tangential to the radial component of the 
proper motion dispersion profile. The solid line represents the model profile, open circles the observational 
data from van Leeuwen et al. (2000), and the black dots the data from Anderson & van der Marel (2010). 




tial relaxation, and brings out the interplay be- 
tween anisotropy and rotation. In fact, by ana- 
lyzing the anisotropy profile derived from the 
observed proper motions (Fig. 3), we note a 
degree of radial anisotropy higher than pre- 
dicted. In fact, it had already been noted that 
radially-biased anisotropic models appear to 
perform better in this cluster (for example, see 
the application of the f^''' models by Zocchi et 
al. 2012 and rotating Wilson 1975 models by 
Sollima et al. 2009). Yet, our model incorpo- 
rates the feature of tangential anisotropy ob- 
served in the outer parts (considered to be a 
natural result of the dynamical evolution of a 
stellar system within an external tidal field) and 
previously exhibited only by orbit-based dy- 
namical studies (van de Ven et al. 2006). 

5. Conclusions 

We have presented a self-consistent global in- 
terpretation of the structure and dynamics of 
0) Cen and 47 Tuc, with particular attention to 
internal rotation, relaxation conditions, and de- 
viations from spherical symmetry. High prior- 
ity has been given to the interpretation of the 
available three-dimensional kinematic data. 

Internal rotation has been shown to play an 
important role in determining the structure of 
the clusters. The difl'erent degrees of success of 
our models have been argued to correspond to 
the different relaxation conditions of these sys- 
tems: the discrepancies noted for co Cen are ex- 
plained by the interplay between rotation and a 
high degree of radial anisotropy. 



The use of proper motions in a dynami- 
cal analysis has been found to be of primary 
importance. By matching the observed proper 
motion profiles with the profiles predicted by 
our models, we have obtained dynamical esti- 
mates of the distances to the two clusters. 
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